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Abstract 

An analytical inactivation cross section formula based on the analytical 
formula for calculating the radial dose distribution of Awad et al. (Awad et 
al., 2018) was proposed. The formula is a multi-hit model based on the track 
structure theory of Katz that was not updated since Katz and co-workers in 
1970s-1980s. The formula was solved numerically to calculation the 
inactivation cross section, o and action cross section for double strand break, 
g D sb for the two strains of the vegetative E. coli cells (B S -i and B/r) 
bombarded with heavy ions of O, Ne, Ar, Kr, Xe, Au, Pb and U of energy 
between 1.1 to 19.4 MeV/n. The response of the cell and hence the 
corresponding inactivation cross section is affected by m, the number of 
targets that must be inactivated per cell; D 0 , the characteristic x or y-ray 
dose; a 0 , the target radius and R max , the maximum range travelled by the 
liberated 8-ray in the medium. The variations in the inactivation cross 
section in terms of these parameters were studied. The model predicts the 
available measured inactivation cross sections data for E. coli cells (B S -i and 
B/r), g as well as the double strand break, g DS b ones with good accuracy. It 
was found that as the ion charge increases as the inactivation probability 
increases and the corresponding inactivation cross section of the ion 
increases. 
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Introduction 


Effect of radiation on DNA and biological tissues is old/new subject 
and is always gain interest in the past, present and future. In radiation 
biology, the amorphous track model was introduced in 1960s of the last 
century (Butts and Katz, 1967; Kobetich and Katz, 1968; Katz et al., 1971) 
and its reduction model, the local effect model (LEM) by Scholz and Kraft 
(Scholz and Kraft, 1994; Scholz and Kraft, 1996; Scholz et al., 1997) were 
used to calculate the biological effect of charged particles. The two 
biophysical models are widely accepted and have been used to calculate the 
effect of radiation on cells and inactivation cross section of ions to enzymes 
and viruses in the last few decades. The two approaches are almost similar 
and the difference between them mostly lies in the concept of local dose 
(Scholz and Kraft) and average dose (Katz) and the definition of the target 
size. In LEM cell inactivation is assumed to be due to production of lethal 
events of certain number in a specific sensitive location (the cell nucleus). 
No interactions of sub lethal damages over large distances of the order of 
micrometers are required to produce lethal events. The present calculations, 
however, is based on the ion kill-gamma kill model of Katz. In which, the 
effects produced by secondary electrons from y-rays and those from 
secondary electrons from heavy ion (8-rays) are comparable at the same 
dose. Two radiation actions were incorporated, y-kill and ion-kill and this 
makes it is possible to explain the difference between the effect of low and 
high-LET radiation. 

Modeling and predicting the radial dose D(r,R) as a function of the 
radial distance, r from the ion pass and 5-ray range, R in a given medium is 
essential for calculating the ion inactivation cross section. The radial dose 
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distribution due to ion interaction with a given medium, D(r,R) was first 
introduced by Butts and Katz in the 1960s (Butts and Katz, 1967; Kobetich 
and Katz, 1968) using a simple empirical electron range-energy relationship. 
This pioneer approach by Katz and co-workers led to formulate many 
biophysical models. It is worth noting that calculating D(r,R) models had 
many modifications and improvement during the last few decades (Katz, 
1978; 1983; 1985; Zhang et al., 1985; Waligorski et al., 1986; Kiefer and 
Straaten, 1986; Cucinotta et al., 1995; Zhang et al., 1994; Chan and Kellerer, 
1997; Cucinotta et al., 1997; Cucinotta et al., 1999, Korcyl, 2012; 
Waligorski et al., 2015 (a,b) and Awad et al., 2018). For more information, 
please refer to these references and references therein. Radial dose, D(r,R) is 
utilized for radiation transport software (Bernal et al., 2015), developing a 
treatment planning system and estimating the cell survival rate in the 
treatment planning system for heavy particle cancer therapy (Waligorski et 
al., 2015-b; Nikjoo et al., 2016) and predict the cell surviving rate and the 
relative biological effectiveness (RBE) of the ions (Waligorski et al., 2015- 
a). 

Measuring the inactivation cross section, o induced by heavy ions on 
viruses, bacteria and DNA were carried out by many authors in the 
literature. For SV-40 virus in EO buffer, o was carried out experimentally as 
well as theoretically (Katz and Wesely; 1991). Inactivation and mutagenic 
effect in spores of Bacillus subtilis ( B. Subtilis ) irradiated with heavy ions 
were investigated by Baltschukat and Florneck (Baltschukat and Horneck, 
1991). Dose dependent decline in surviving fraction in Chines hamster V79 
cells exposed to high LET ions has been studied by Pathak et al. (Pathak et 
al., 2007). Experimental and theoretical cross sections for E. coli mutants B, 
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B/r, and B s _i after irradiated by Heavy-Ion Irradiation were studied (Katz 
and Zachariah, 1993). Inactivation cross section of ions for dry enzymes and 
viruses were calculated using an empirical formula (Liu et al., 1994 and Liu 
et al., 1996). Heavy ions induced inactivation and mutagenic effect in 
microorganisms such as yeast cells which provide astronauts needs in active 
space exploration was studied (Wang el al., 2010). Quantitative analysis of 
radio-induced DNA damage exposed to protons at Bragg-peak energy was 
given by Souici et al., (Souici et al., 2017). 

The aim of the present work is to calculate the inactivation cross 
section, o using a new analytical formula. This formula is based on the 
analytical approach proposed by Awad et al. (Awad et al., 2018) for 
calculating the radial dose distribution D(r,R). The formula is a multi-hit 
model based on the track structure theory of Katz and is considered as an 
alternative to Monte Carlo code for calculating the radial dose distribution in 
a given medium due to the passing charged ion. The calculations were 
carried out for two strains of Bacterial cells of E.Coli. Estimating the 
inactivation cross sections, o for the two strains (B s _! and B/r) E. coli cells 
bombarded with heavy ions of O, Ne, Ar, Kr, Xe, Au, Pb and U of energy 
between 1.1 to 19.4 MeV/n were carried out. In addition to calculate the 
induction double strand break cross section, o DS b for those two strains. The 
different parameters that affect a of m, D 0 , a 0 and R max were studied. The 
calculated inactivation cross sections were as compared to the measured 
ones. The inactivation probability for each ion was studied. 
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Experimental data 


The present calculations were compared to data that obtained from a well- 
designed experiment in which different ions of different charges with almost 
similar energy/n were utilized to irradiate the samples. Samples of two 
strains (B s _! and B/r) of the vegetative cells of Escherichia coli (E. coli) of 
almost similar DNA structure were used (Schafer et al., 1994). They are 
considered as a simple test system to measure the induction of DNA double 
strand breaks (DSB) and their correlation with cell survival dependent on the 
radiation quality. The two E. coli strains were bombarded with heavy ions of 
Ne, Ar, Kr, Au, Pb and U. The inactivation cross section and double strand 
break cross section were measured. Another inactivation cross section data 
for B S -i and B/r E. coli bombarded with O, Ne and Xe were obtained 
(Schafer et al., 1987) were added. Table 2 compile the different ions of 
energy between 1.1 to 19.4 MeV/n quoted in the present work and there 
references. Schafer et al 1994 data will be superimposed in the present work 
as (Exp), while Schafer et al., 1987 data will be quoted as (Exp report). 


Methodology 

Inactivation (loss of colony forming ability) of enzymes, bacteria and 
viruses by energetic heavy ions gained large interest and still, therefore, 
many approaches were suggested for calculating the inactivation cross 
section. Calculating inactivation probability, P and a is mainly depends on 

D(v R 1 

the —— ratio. For each D(r,R) formula there is inactivation cross section 

D o 

formula. Therefore, there are Katz (Butts and Katz; 1867), Zhang (Zhang et 
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al., 1985); Waligorski (Waligorski et al, 1987); Cucinotta (Cucinotta et al., 
1997) as well as Waligorski (Korcyl, 2012; Waligorski et al., 2015) 
approaches. The present calculations (will call it Calc) work will be 
compared with Katz and Zhang calculations (Zhang et al., 1985) and will 
call it (Katz and Zhang). 

The multi-hit inactivation cross section is a statistical model 
introduced by Katz (Katz et al., 1971; Katz et al., 1996). In that model it was 
assumed that the effects produced by secondary electrons from y-rays and 
those from secondary 8-electrons from heavy ions are comparable at the 
same dose. Therefore, the probability for producing effects in a macroscopic 
volume by a given dose of y-rays in the target was used to estimate the 
probability of producing effects in that target by a charged ion as follows: 

r D(r,R) \ 

p = l - e ( D 0 ) (!) 

where, D 0 is the y-ray dose for 37% survival of the cell and D(r,R) is the 
average radial dose over the target volume. D(r,R) is assumed to be the 
energy deposited per unit mass in a short cylinder of radius a 0 , whose axis is 
parallel to the path of the ion. A semi-empirical analytical model based on 
electronic radiation damage is introduced by (Awad et al., 2018) for 
calculating the radial dose distribution D(r,R), In this approach, the 
empirical electron range-energy formula of Tabata et al. (Tabata et al., 1972) 
was considered where the electron range, R in cm is given as: 

R( r) = a 1 (^n(l + o 2 ^)--^ r ), (2) 

where 
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a ± = 


b x A T 

~pA 2 


q>2 — b 2 Zj, 


a 4 — bjZj’, 


a K — 


a 3 — b 4 b 5 Z T , 
bs 


5 ybq 

Z V rj-1 


where Z T and A r are the atomic and mass number of the target, respectively. 
bi (i=l,2,...,9) are constants independent of absorber material and are given 
in Table 1. Tis the energy of the ejected 5-electrons and its maximum value 
as a function of /3 (the ion velocity in medium relative to velocity of light in 
vacuum) is given by: 

2 mc 2 p 2 


T = 

1 max 


1 -/ 3 2 


(3) 


The second important step is to calculate the electron energy in MeV as a 
function of the electron range in cm. The inverse of the energy-range of Eq. 
(2) was found by Tabata et al. (Tabata, et al., 1972) as follows: 


T(R) = 0.511 c x ^e c i (C2+1+c 4« C5) - l), 


(4) 


where c, (i=l,2...5) are constant for a given target (Z T , A r ) and d t 
(i=l,2,...,9) are another constants independent of absorber material and are 
given in Table 1 as: 



C 2 ~ 


d 2 Zy 3 


Cg — d 4 d 5 Z 


t > 


d 6 _ d 8 

c 4 “ 3*7 ' c 5 “ “ 

Zj rj-i Zi rj-1 

These two steps are important for defining the necessary constants and for 
more details, please refer to Awad et al., 2018. Finally, the radial dose, 
D(r, R ) depends on r and R and its formula was defined as follows: 
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D(r,R ) 


-CZ* 2 e x 4 x 3 (x 2 -Xi) eV c i Xl 
p2nrc 1 (3 2 (e x 4-i)2 


( 5 ) 


Where 


l+c 4 q?-r) c 5 

(6) 

Y _ c 3 c 4 (R-r) c 5 

2 (l+c 4 (7?-r) c 5) 2 

(7) 

r C3 (l + c 4 /? c 5) C 3 C 4 c 5 RC 5 

Y — C 2 1 c 4 C 1 

3 Ci (1 + C 4 R C 5)2 

(8) 

Z 4 =-(c 2 + c ' 3 c ) 

4 C 4 V z 1 + C 4 /? C 5/ 

(9) 


Q7TL 

For water medium, p — 1 and C is a constant (Zhang et al., 1985) equals 
C = ^4- = 8.5— = 1.369 x 10“ 7 — = 1.369 x lO" 14 — (10) 

me 2 mm cm cm 

The inactivation probability thus can be determined by using Eq. (5) for 
D(r, R) and hence, the inactivation cross section, o for a given ion is given 
as follows: 

„ p (D{r,K)\-y m 

o — 271 f max f r 1 - e ^ °o ) dr dR (11) 

^min 1 

where, m represents the number of targets in the cell which must be hit and 
inactivated by radiation to produce an observable effect in such cell. For one 
hit detectors m = 1 i.e. a single hit in a single target suffices to inactivate the 
detector. For biological cells m is frequently 2 or more (Katz et al., 1996). 

An algorithm was constructed to calculate the double integration of o 
as given in Eq. (11) by using Mid-Point numerical Method of integration for 
a given target. For calculating the radial dose, the ion parameters of Z* (the 
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effective charge number of the ion) which was calculated according to 

_2 

Barkas’ formula as Z* — Z( 1 — e _125 ^ z 3 ) and (3 were incorporated in the 
calculation as well as the target parameter of p = 1 where the E. Coli 

density was considered as a water medium. The two integration limits were 
carried out in such a way that the down limit for the second integration is 
R min and the upper limit is R max are the limits for radial dose estimation (Eq. 
5). For the first integration that is given in Eq. 11, the lower limit is taken as 
1 nm and the upper limit is less than or equal R min of the first integration. It 
is worth noting that, the first integration is the target radius a 0 . 

3. Results 

3.1 Parameters affecting inactivation cross section calculations 

The response of the E. coli cells (B S -i and B/r) to the bombarded ions 
was studied through investigating the inactivation cross section of these cells 
using Eq. 11. Equation 11 was solved numerically to calculation the 
inactivation cross section, o and induction double strand break cross section, 
a dsb for the two strains of the vegetative cell of E. coli cells (B s _! and B/r) 
bombarded with heavy ions of O, Ne, Ar, Kr, Xe, Au, Pb and U of energy 
between 1.1 to 19.4 MeV/n. The data for each strain will be demonstrated 
and discussed separately later on. 

Inactivation cross section depends on four parameters (m, D 0 , a 0 and 
Rmax)- The number of targets that must be hit by radiation to cause an 
observable effect (inactivate) per cell is m. D 0 is the characteristic dose and 
it is the y-ray or x-ray dose for 37% of cell survival, a 0 is the target radius 
and finally, R max is the maximum 5-electrons range corresponding to the 
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maximum 8-ray energy transferred due to the energetic charged ion (Eq. 3) 
and it is the upper limit for the second integration in Eq. 11. In this part, the 
effect of these parameters on o will be discussed and inactivating of B/r 
E.coli cells by Argon ion will be presented as an example. The variation in 
o with the ion energy at different m values of m=l ,2,3,4 and 5 and at D 0 =61 
Gy and target radius, a 0 =0.2 pm is shown in Fig. 1. One can observe that, the 
inactive cross section decreases with m and m=2 gives the best fit in this 
case. It was found that o is very sensitive to m and m=2 give the best fit for 
g for the two E. coli cells (B and B/r) and m=3 give the best fit for the 
double strand break action cross section, g DS b for E. coli cell DNA. 

Inactivation cross section as a function the Ar energy for different 
target radius of a 0 =0.1, 0.2, 0.3 and 0.5 pm and at m=2 and D 0 =61 Gy is 
illustrated in Fig.2. It can be seen that, at low energy (less than ~8 Mev/n) 
the target radius is slightly effective while for energy greater than 8 MeV/n 
is a 0 is not changing o significantly. It was found that target radius that gives 
best fit to the experimental data is at a 0 = 0.2 pm (200 nm). Therefore, all 
targets radii will be fixed at 0.2 pm. 

Regarding the characteristic dose, D 0 the inactivation cross section for 
B/r E. coli cell was studied at different Ar energy and at different D 0 values, 
see Fig. 3. One can observe that D 0 values are significantly changed the 
obtained inactivation cross section and o is very sensitive to D 0 . In opposite 
to previous study, the present work was fitting the inactivation cross section 
using different D 0 values for each ion. The previous studies, were fixing D 0 
value for all ions. Table 3 is gathering the characteristic dose, D 0 for the 
different ions in comparison with the literature D 0 values for both vegetable 
E. coli cells. 
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Table 4 as well as Fig.4 demonstrates the effect of R max , the upper 
integration limit on inactivating B/r E. coli cells at different Ar energy. 
Using the full 8-ray range (R max ) that was deduced using Eq. 3 and Eq. 2 
gives progressive cross section values. It was found that at energy between 2 
to 10 MeV/n the deduced R max using Eq. 3 and Eq. 2 fits most of the studied 
data, however, a problem exist in the region above 10 MeV/n where the 
calculated o cannot fit the experimental data at such high energy regions. 
Therefore, many efforts and by using trial and error method, the present 
work suggested R max values that fit most of the studied experimental data 
and is given in Table 4. 

3.2 B S -i E. coli inactivation cross section 

In the energy range 2-20 MeV/n, the present model calculated the 
inactivation cross section for B s .i E. coli cells bombarded with O and Ne 
ions using the fitting parameters m=2, a 0 =0.2 pm and D 0 = 12.5 Gy (see 
Table 3). o (pm 2 ) as a function of O and Ne energy was plotted in Fig. 5. 
Good agreement between the present calculations and the experimental data 
for O (O, Exp.) and Ne (Ne, Exp.) data and it is closer to Zhang and Katz 
calculations for O as well. Similarly, the experimental data for B S -i E. coli 
cells strike by Ar and Kr at different energies are well predicted by the 
present model, Fig.6. Table 3 compiles the implemented model parameters. 
The present calculations as well as the experimental data are in agreement 
with calculated data by Katz and Zhang for both ions. 

Fig. 7 and Fig. 8 present a comparison between experimental 
inactivation cross section for B s .i E. coli cells bombarded by Pb, Au, Xe and 
U ions. The fitting parameters for these ions are compiled in Table 3. One 
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can observe that as the ion charge increases as the cross section of 
inactivation increases. Reasonable matching between the present 
calculations and the experimental data at both low and high energy part in 
agreement with Katz and Zhang calculations as well. 

3.3 B/r E. coli inactivation cross section 

The present model is able to predict the inactivation cross section for 
B/r E. coli cells bombarded by heavy ions quit well in comparison with both 
experimental and Katz and Zhang calculations. It must emphasis that 
inactivation cross section for B/r E. coli cells was carried out using target 
radius, a 0 =0.2 pm and m=2. Comparison of the predicted inactivation cross 
section for B/r E. coli cells produced by the present approach and the 
experimental inactivation of B/r E. coli cells bombarded by O and Ne is 
illustrated in Fig. 9. The corresponding fitting parameters for this strain are 
given in Table 3. Inspecting Table 3 shows that B/r E. coli cells are more 
resistant to radiation than the B s _ 7 E. coli cells. It was found that D 0 for B/r 
E. coli is always larger than D 0 for B S -i E. coli cells at a given ion and at a 
given energy. 

Inactivation of B/r E. coli cells by Argon, Krypton, Xenon, Lead and 
gold ions as a function of energy/n are shown in Fig. 10 and Fig. 11. General 
satisfactory agreement has been found among the present calculations, 
experimental data and Katz and Zhang calculations especially at the low 
energy part. However, some deviations appeared between calculated and 
experimental data at high energy part. 
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3.4 Action cross section for DSB 


Interaction of radiation with DNA causing double strand break is the 
major lesions in the genomes, therefore, the action cross section for DNA 
samples, g DS b for E. coli cells (B s _i and B/r) was also calculated using Eq. 
11. The beauty of the present approach is its ability to predict the action 
cross section for double strand break, g DS b for E. coli cells (B s _! and B/r) 
DNA using the same fitting parameters (Table 3 and 5) for a given strain or 
cell. The only difference is that DNA needs three hits (m=3) while the cell 
itself is inactivated by two hits only ( m=2 ). DNA samples were bombarded 
with different ions at different energies. The different fitting parameters for 
each ion are compiled in Table 5. It was found that the experimental data 
was well predicted using the target radius, a 0 =0.2 pm and m=3 and D 0 values 
are different for different ions. In this table, the experimental action cross 
section for double strand break, g DS b for E. coli cells (B s _i and B/r) DNA are 
compared to the calculated ones. Global agreement between experimental 
and calculated data was found especially at the low energy ions but at higher 
energy still some disagreement between calculations and experimental data. 

3.5 Ion’s inactivation probability 

The inactivation probabilities for B s .i E. coli, as an example, 
bombarded with a series of different ions were studied. Inactivation 
probability (Eq. 1) of some charged particles of different atomic number 
(Ne, Ar, Kr, Au, Pb and U) having similar energy per nucleon was 
calculated as a function of the radial distance, r normal to the ion pass. 
Inspecting Fig. 12, one can observe that for light ions i.e. Neon can 
inactivate the bacteria -100% if it hit it within - 1 nm from the ion path 
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while Argon ion at ~ 2 nm, Krypton at ~ 5 nm. Heavy elements like lead and 
U can inactivate the bacteria a 100 % if it hit it within ~ 15 nm and Au is 
within 20 nm from the ion projectile. 


4. Discussion 

Inactivation cross section of two E. coli strains (B s _! and B/r) by 
heavy elements was calculated by integration the inactivation cross section, 
Eq. 11 numerically by the mid-point method. Four parameters are affecting 
the inactivation cross section values i.e. o (a 0 , m, D 0 , R max ). It was found 
that inactivation cross section is slightly depended on the target radius a 0 and 
a 0 equal 200 nm or 0.2 pm fits the experimental data for all ions. Therefore, 
a 0 =0.2 pm kept constant for o and o D sb through-out the present calculations 
for B s _! and B/r cells. This target radius is comparable to the radius of the 
cell where E. coli bacterium is about 1-3 pm long and 1 pm width and about 
0.25 pm in diameter. The target radius, a 0 equal 0.2 pm is closer to previous 
studies, see Table 3. 

The number of targets that must be hit by radiation to cause an 
observable effect (inactivate) per cell was found equal 2 (m=2) by the 
present calculation for E. coli cells (B s _i and B/r). In the meantime, in 
calculating the action cross section for DNA for the same E. coli (B S -i and 
B/r) cells, g DS b it was found that m=3 give the best fit to the experimental 
data. Therefore, E. coli cells (B s _! and B/r) are not 1-hit detector (m=l) as 
Katz calculated (Katz, 1996), instead it is a multi-hit detector as expected 
before by Caucinotta et al., 1996 and Katz 1996 for biological systems. This 
implies that two electrons passing through the bacterium (B s _i and B/r) and 
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three electrons for their DNA are capable of inactivating E. coli cells. E. coli 
cells are multi-hit detector is more reasonable compared to the notion that 
it's one-hit system (a single electron passing through the bacterium is 
capable of inactivating the cell). For biological cells m is frequently 2 or 
more (Katz et al., 1993 and Katz et al., 1996). The present finding of m=2 is 
more reasonable in this regard and suggest the significant of the direct 
interaction between the incoming ions and cells rather than the indirect 
interaction. 

The characteristic dose, D 0 is significantly affecting the calculated 
inactivation cross section; a. Katz and co-workers were assuming fixed a 0 
and D 0 values but in the present work a 0 was fixed at 0.2 pm while D 0 was 
varying with different ions. Table 3 and compiles the D 0 values for E. coli 
strains and its DNA for the different bombarded ions. One can observe that 
D 0 increases when the atomic number increases as stated in Tables 3 and 5. 
It was found that the results of the present calculations depend on D 0 than 
the target radius, a 0 in agreement with what was observed by Waligorski et 
al., 1987. Fitting the experimental data, suggested the change in D 0 gives 
better prediction to the experimental data rather than fixing it as previous 
investigator do. B/r E. coli is more resistant to radiation compared to B S -i 
strain where D 0 for it is always large as illustrated in the tables. 

Inspecting the figures as well as the tables, one may observe that the 
calculated inactivation cross section for heavy ions is not 100% in agreement 
with the experimental data. One of the probable reasons for this 
disagreement could be the mean dose distribution over the target radius 
when particle are passing the target area (Schafer et al., 1994 and Katz et al., 
1996). This assumption is a fundamental one suggested by Katz in his multi- 
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hit model. In comparison with experimental data, the inactivation cross 
section has the difficulty due to determining the real maximum range, R max 
travelled by 5-rays (inaccurate determination for R max ). R ma x is the upper 
integration limit and it is significant in both calculating the radial dose 
(Eq.5) and calculating o (Eq. 11). Therefore, calculating the maximum 5-ray 
range by Eq. 3 and Eq. 2 give good fit to the experimental data up to 10 
MeV/n. But larger than 10-12 MeV/n, there are some disagreement between 
the experimental and the calculated data. Therefore, the trial and error 
method was applied and the suggested R max that improved the model 
prediction and reasonable agreement with experimental data was achieved. 
This new R max can be accepted in the view that electrons are not moving in 
straight line and it suffers from multiple scattering in its path. One important 
reason for this slightly disagreement at higher energy parts may be the 
scattering in the experimental data itself. In addition to this, Takahashi et al., 
1983 in his attempt to explain the disagreement between experimental and 
calculated data (Takahashi et al., 1983) suggested that closer to the ion path 

o 

(10 A) most atoms or molecules are ionized in such a way that such densely 
excited or ionized species may react with each other and cause thermal spike 
effects. Such thermal effect may cause derangement of biological systems 
and produces boost growing of inactivation cross sections within this narrow 
region closer to the ion. Therefore, the mean dose distribution over the 
whole target radius assumption may not exactly hold. 
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5. Conclusion 


The present approach is a promising one for calculating the 
inactivation cross section for bacteria and viruses where good global 
agreement between the present calculations and the experimental data was 
obtained. 

Inactivation cross section calculation and cell death probability due to 
radiation is the target of this work and it was found that ions with higher 
atomic numbers produce larger inactivation probability and inactivation 
cross section to the E. coli cells (B s _! and B/r) and its DNA compared to the 
light ions. 

E. coli cells are a multi-hit detector where m=2 for the two E. coli 
strains and m=3 for its DNA rather than a 1-hit detector. 

Calculating the inactivation cross section for B/r E. coli cells shows 
that B/r E. coli is more resistant or less sensitive to radiation than the B s _, E. 
coli cells where the characteristic dose D 0 for B/r E. coli is always larger 
than D 0 for B s .i E. coli cells at a given ion. 

Calculated inactivation cross section is slightly depending on the 
target radius, a 0 . In the meantime, it is strongly depend on the chosen R max , 
the travelled 8-ray range in the medium. 

The present model predict inactivation cross section for lighter ions 
(O, Ne, Ar and Ke) very well but for heavier ions (Xe, Au, Pb and U) it fit 
the experimental data very well at lower energies (less than 10 MeV/n) but it 
found some difficulties to predict the data at higher energies. The scattering 
in the experimental data helped to increase this observation. It may be 
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attributed to inaccurate R max determination and the condensed thermal spike 
region produced by the heavy ions. 

In progress, calculating the inactivation cross section for the more 
complicated mammalian cells; animal and culture cell using the current 
approach. Estimating the relative biological effectiveness, RBE for the 
different ions as well as using this method for developing treatment 
radiotherapy plans for some ions like carbon is under investigations. 
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Figure Captions 

Fig.l Calculated inactivation cross section at different m values for 

B/r E. coli cells bombarded by Argon at different energies in 
comparison with experimental data. 

Fig.2 Calculated inactivation cross section at different target radii, a 0 

for B/r E. coli cells bombarded by Argon at different energies 
in comparison with experimental data. 

Fig.3 Calculated inactivation cross section at different characteristic 

dose, D 0 for B/r E. coli cells bombarded by Argon at different 
energies in comparison with experimental data. 

Fig.4 5-ray maximum range, R max effect on the calculated 

inactivation cross section for B/r E. coli cells bombarded by 
Argon at different energies in comparison with experimental 
data. 

Fig.5 Inactivation cross section, o (pm 2 ) of B s _i E. coli cells by O and 

Ne ions at energy 2-20 MeV/n in comparison with experimental 
data as well as Katz and Zhang calculations. 

Fig.6 Inactivation cross section, o (pm 2 ) of B s _! E. coli cells by Ar 

and Kr ions at energy 2-20 MeV/n in comparison with 
experimental data as well as Katz and Zhang calculations. 
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Fig.7 Inactivation cross section, a (pm 2 ) of B s _i E. coli cells by Pb 

and Au ions at energy 2-20 MeV/n in comparison with 
experimental data as well as Katz and Zhang calculations. 

Fig.8 Inactivation cross section, o (pm 2 ) of B s _] E. coli cells by Xe 

and U ions at energy 2-20 MeV/n in comparison with 
experimental data as well as Katz and Zhang calculations. 

Fig.9 Inactivation cross section, o (pm 2 ) of B/r E. coli cells by O and 

Ne ions at different energies in comparison with experimental 
data as well as Katz and Zhang calculations. 

Fig. 10 Inactivation cross section, o (pm 2 ) of B/r E. coli cells by Ar and 
Kr ions at different energies in comparison with experimental 
data as well as Katz and Zhang calculations. 

Fig.l 1 Inactivation cross section, o (pm 2 ) of B/r E. coli cells by Xe, Pb 
and Au ions at different energies in comparison with 
experimental data as well as Katz and Zhang calculations. 

Fig. 12 Inactivation probability of B s _! E. coli cells bombarded by Ne 
(10.48 MeV/n), Ar (10.49 MeV/n), Kr (11.22 MeV/n), Au 
(10.33 MeV/n), Pb (12.75 MeV/n) and U (13.52 MeV/n) ions 
as a function of the radial distance, r from the ion trajectory. 
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Table Captions 

Table 1 Values of the constants b t and used for electron range-energy 
equation, Eq. (1) and its inverse relation of energy-range 
equation, Eq. (2). 

Table 2 The experimental data are mainly from Schafer et al., 1994 

(Exp.) and Schafer et al., 1987 (Exp. report) that will be marked 

* 

as . 

Table 3 The implemented characteristic dose, D 0 for the different ions 
in comparison with literature value for B S -i and B/r E. coli cells. 

Table 4 The upper integration limit, R max where the lower integration 
limit R min =1 nm. 

Table 5 Experimental inactivating Double Strand Break cross section 
(o D sb) data of E. Coli cells (B S -i and B/r) bombarded by 
different ions in comparison the calculated o DS b ones using 
m=3. 
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Table 1 


i 

bi 

di 

1 

0.2335 ± 0.0091 

(2.98 ± 0.3)xl0 3 

2 

1.209 ± 0.015 

6.14 ± 0.29 

3 

(1.78 ± 0.36 )xl0 -4 

1.026 ± 0.02 

4 

0.9891 ± 0.001 

(2.57 ± 0.12)xl0 2 

5 

(3.01 ± 0.35)xl0' 4 

0.34 ± 0.19 

6 

1.468 ± 0.09 

(1.47 + 0.19)xl0 3 

7 

(1.18 ± 0.097)xl0 2 

0.692 ± 0.039 

8 

1.232 ± 0.067 

0.905 ± 0.031 

9 

0.109 ± 0.017 

0.1874 ± 0.0086 
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Table 2 



Bj 





Table 3 


Literature 


Present 

D 0 (Gy) 

Do (Gy) 

References 

a 0 (|um) 

Bs-i 

B/r 

a 0 =0.2 

jam 

Bs-i 

B/r 

Z 

o Odsb 

o Odsb 

Takahashi et al., 
1986 


12.6 

36.5 

8 

12.5 

20 

Schafer et al., 1994 

0.2/0.4 

15.4 

47.6 

10 

20 

29 

Zhang and Katz, 
1995 

0.5 

13.8 

44.6 

18 

42 

61 

Katz et al., 1996 

0.1 

12.6 

40 

36 

70 

80 

Liu et al., 1996 


12.6 

47.6 

54 

77 

110 

Korcyl, 2012 

0.1 

12.6 


79 

110 

185 

• 




82 

120 

170 





92 

100 
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Table 4 


E/n 

(MeV/n) 

^max (Will) 

(Eq. 3&Eq.2) 

present work 

2 

660 

660 

4 

1930 

1930 

6 

3800 

3800 

8 

6240 

6240 

10 

9240 

9240 

12 

12800 

11500 

14 

16800 

12800 

16 

21400 

12800 

18 

26400 

12800 

20 

31600 

12800 


31 




Table 5 


z 

Energy 

(MeV/n) 

Do (Gy) 

DSB 

Bs-i 

B/r 

Bs-i 

B/r 

^Exp 

(bm 2 ) 

^Calc 

(|um 2 ) 

^Exp 

(lunr) 


10 

14.09 

20 

29 

0.41+0.04 

0.29 

0.26+0.08 

0.13 


11.77 

20 


0.36+0.06 

0.42 




10.48 

20 

29 

0.44+0.08 

0.41 

0.10+0.07 

0.17 


4.43 

20 

29 

0.41+0.14 

0.35 

0.31+0.08 

0.16 

18 

14.36 

42 

61 

0.66+0.06 

0.57 

0.36+0.06 



13.59 

42 

61 

0.70+0.11 

0.65 

0.40+0.11 



10.49 

42 


0.81+0.12 

0.73 




7.6 

42 

61 

0.51+0.20 

0.65 

0.35+0.27 



4.7 

42 

61 

0.63+0.13 

0.52 

0.73+0.15 


36 

13.22 

70 


0.72+0.18 

2.12 




11.22 

70 

80 

1.58+0.16 

2.35 

1.78+0.08 



6.24 

70 

80 

1.19+0.09 

1.48 

1.29+0.09 


79 

10.33 

110 

185 

3.06+0.46 

9.35 

3.49+0.33 



5.1 

110 

185 

2.81+0.67 

3.85 

2.54+0.47 


82 

12.75 

120 

170 

4.10+0.55 

6.20 

3.25+0.81 



4.28 

120 

170 

3.20+0.45 

1.90 

2.22+0.43 


92 

13.52 

100 


4.59+0.97 

10.5 
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